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CONVECTION OF A PATTERN OF VOETICITI THROUGH A SHOCK WAVE 

By H. S, Ribner 


SUMMAEI 

An arbitrary weak spatial distribution of vorticity can be repre- 
sented in terms of plane sinusoidal shear waves of all orientations and 
wave lengths (Fourier integral). The analysis treats the passage of a 
single representative weak shear wave through a plane shock and shows 
refraction and modification of the shear wave with simultaneous gener- 
ation of an acoustically intense sound wave. Applications to turbulence 
and to noise in supersonic wind tunnels are indicated. 


INTRODUCTION 

Turbulence such as the residual small eddying motion in a wind- 
tunnel stream will gradually decay as it is carried along. The decay 
process has been the subject of much study in the face of formidable 
difficulties. The random character of the motions has been success- 
fully handled by the methods of statistics; even with these methods, 
however, the nonlinearity of the eq^uations governing the intennixing 
processes has severely limited the progress attainable without 
simplifying assumptions. 

On the other hand, for relatively sudden changes in turbulence, 
such as occur when it passes through a wire-mesh damping screen ,the 
decay may be negligible and the changes may follow linear laws. The 
linearity is assured if the turbulence constitutes a sufficiently 
small perturbation of the main stream. Recently it has been found 
that the problem of such linear changes could be solved completely by 
a specialized adaptation of the spectrum concept of the statistical 
theory of turbulence. 

Several of these linear processes have been treated in thld 
manner: the damping-screen problem (reference l), the passage of 

turbulence through a sudden wind-tunnel contraction (reference 2), 
and the passage of turbulence through a series of screens followed 
by a sudden contraction (unpublished investigation of M. Tucker) . A 
basic techniq.ue for such problems has been evolved in these papers. 
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The present paper is motivated hy another problem of the same linear 
character, namely, the convection of veak turbulence through a shock wave. 
Among other circumstances, this problem arises in the interpretation of 
measurements with a hot-wire anemometer in a supersonic stream, because 
a detached bow wave stands ahead of the wire.^ Such a curved shock is 
not attractive for theoretical analysis, but it is not difficult to 
replace it with an extended plane shock by use of auxiliary means; 
attention can thus be limited to the convection of turbulence through 
a plane shock. 


The conceptual basis for the treatment of these linear problems is 
as follows: An arbitrary weak spatial distribution of vorticity - and 

hence a weak turbulent velocity field - can be represented as a super- 
position or spectrum of plane sinusoidal shear waves distributed among 
all orientations and wave lengths. This is a physical interpretation 
of the mathematical formulation as a Fourier integral;^ the individual 
shear waves may be identified as Fourier or spectrum components. When 
the turbulence wave pattern is convected through a screen or through 
a shock wave, the individual waves are altered without mutual inter- 
ference if the waves are suitably weeik. Thus the modified field down- 
stream of the screen or shock can be obtained, in principle, by super- 
position of the modified Individual waves. In practice the description 
of the detailed spatial distribution of velocity, either initially or 
finally, is hopeless; the initial wave distribution is known only 
statistically (e.g., the phase angles are unknown), and statistical 


^A simple interpretation for all but very small eddies comparable 
with the scale of the bow wave is, however, available in the work of 
Kovasznay (reference 3). 

^The velocity field so represented may be either rotational or 
irrotational within the specified region, even though the building 
blocks," the shear waves, are rotational. In case an irrotational field 
is represented, the vorticity of these shear waves, but not the veloc- 
ity, mutually cancels within the specified region (which may be mul- 
tiply connected), leavlsig a distribution of vorticity in the external 
space. The irrotational flow may be regarded as induced by this 
external vorticity. 

These remarks all refer to a velocity field satisfying the incom- 
pressible continuity equation: a small- perturbation field of vorticity 

in fluid at rest, or convected by a main stream, will fulfill this con- 
dition. 
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changes only can be calculated. In either case the analysis of the 
behavior of a representative single wave constitutes a prerequisite to 
the detenaination of the changes in the weak turbulent field. 

In the present paper such an analysis is carried out for a single 
shear wave, of arbitrary inclination, convected through a plane shock. 
There remains the task of calculating therefrom the changes in the 
statistical properties of a weak turbulent field convected through 
a shock. ^ 

This single-wave problem is also treated in a current investigation 
by F. K, Moore (unpublished). The analyses bear little resemblance: 

In that work a reference frame is used in which the flow is unsteady, 
whereas herein a frame is used in which the flow is steady, Sound 
waves are likewise treated in the work cited. 

The outline of the present analysis is as follows: The problem is 

posed as the calculation of the flow field behind a plane normal shock 
wave due to the convection throu^ the shock of an inclined plane sinu- 
soidal shear wave; the shear wave is specified to be weak to ensure 
small perturbations to the mean flow. This problem, for which the 
flow is unsteady in time, is converted into an equivalent steady-flow 
problem by transformation to a moving frame of reference. In this 
frame the normal shock is replaced by an equivalent oblique shock. 

The analysis is now formulated as a boundary-value problem for the 
flow in the region downstream of the shock: The governing partial 

differential equation for this small-perturbation rotational flow is 
derived (extension of Sears* work, reference 4); boundary conditions 
on the velocity components Just behind the shock are obtained from the 
oblique-shock relations; and finally the rotation term in the governing 
equation is evaluated in terms of gradients of entropy and total 
enthalpy, with use of the entropy changes across the shock. The 
initially unknown perturbation of the fonm of the shock wave is taken 
into account in the boundary conditions and rotation term by assuming 
it to be sinusoidal with initially undetermined amplitude eind phase. 

The velocity W (all symbols are defined in appendix A) down- 
stream of this equivalent oblique shock may be either subsonic or 
supersonic depending on the inclination of the initial plane sheeir 
wave; separate solutions of the boundary- value problem are worked out 
for the two markedly different cases. The horizontal shear wave - 
which is a simple special case for subsonic W - is given a separate 
treatment . 


^Procedures have been developed in references 1 and 2; an exten- 
sion will be required if the noise field generated by the interaction 
is to be treated. 
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The analysis of the velocity field downstream of the shock is fol- 
lowed hy an account of the associated pressure and density fields there 
and of the distortion of the initially plane shock. Finally, the 
acoustic level of the (fluctuating) pressure field is worked out in 
approximate fashion for an example applied to a supersonic wind 
tunnel; A particular initial intensity of turbulence is assumed and 
considered as being concentrated in a single shear wave rather than 
distributed throughout a continuous spectrum. The calculation amounts 
to an estimate of the noise level generated by the passage of a 
specified level of turbulence through a shock wave. 

This investigation was conducted at the NACA Lewis laboratory. 


FOEMUIATION OF BOUKDAEI-VALDE PEOBLEM 

The unsteady-flow problem . - The Inclined plane sinusoidal shear 
wave is shown schematically in figure 1. The flow is viewed in a plane 
perpendicular to the shock and to the wave fronts. The wave is supposed 
to be convected downstream by the main stream with velocity so that 

it passes through the normal shock. The passage through the shock is 
evidently an unsteady process, since the intercepts of the Inclined 
lines (the nodes of the sine wave) move downward along the shock 
front; it will be shown that a disturbance ripple moves along the 
shock with the same speed V. 

In the general' case of a plane obliq.ue sinusoidal shear wave there 
will also exist a perturbation velocity component normal to the plane 
of the figure. Now the ripples in the shock front will be two dimen- 
sional, and the shock with the ripples will still be everywhere perpen- 
dicular to the plane of the figure. Thus, the normal velocity component 
will be parallel everywhere to the shock aind will be unaffected as the 
shear wave passes through; the component will have no other effect. Its 
invariance established, this normal velocity component will be omitted 
from the analysis. 

The equivalent steady-flow problem . - If an observer moves down- 
ward along the shock with a speed Y, relative to him the flow will 
have an apparent upward velocity component V. This scheme of things 
is shown in figure 2, In particular, Y has been chosen so that the 
resultant stream velocity (relative to the moving observer) is alined 
with the velocity in the disturbance wave; that is, Y = tan 0. The 

observer then sees what appears to be a steady sinusoidal shear flow 
passing through an oblique shock wave. Thus, by the proper choice of - 
a system of moving axes the original unsteady-flow problem has been 
converted into an equivalent steady-flow problem. 

Governing partial differential equation for rotational flow . - The 
task of the analysis is to calculate the flow field on the downstream 
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side produced by the passage of the sinusoidal shear flow through the 
equivalent oblique shock. It is to be expected that the shock will be 
perturbed from its mean plane and will, in fact, develop a corrugated 
appearance. Because of these corrugations, vorticity (rotation) will be 
introduced into the downstream flow. This vorticity and all the down- 
stream velocity perturbations will be weak compared with the stream 
velocity because the original disturbance wave has been assumed weak. 
Thus, a small-perturbation, or linearized, treatment of the flow 
field is permissible. 

In reference 4 the governing partial differential equation for 
small-perturbation compreseible rotational flow has been derived for 
isoenergetic flow, that is, for flow of constant stagnation enthalpy. 
However, the shear wave under consideration possesses variable stagna- 
tion enthalpy; that is, pressure, density, and temperature are constant 
upstream of the shock, but the velocity varies. It has been necessary, 
therefore, to obtain a more general governing equation that applies 
when both entropy and stagnation enthalpy are variable. The derivation 
is given in appendix B. 

This governing equation is expressed in terms of coordinates ^ 
and T), 4 being the distance in the main stream direction and q 
the distance perpendicular thereto. The equation reads 

(1 + .a (1) 

where W is the stream velocity in the transformed problem, W is 
the corresponding Mach number, H is the stagnation enthalpy, s is 
the entropy, T is the temperature, 2 is the vorticity, and \|r is a 
perturbation stream function such that 

\|r^ H w = perturbation velocity in ? direction 
-(1 - = w* = perturbation velocity in direction 

(The stream function is defined differently in reference 4, as it 
involves an entropy term. ) 

For application of equation (l) in the present problem reference 
should be made to figure 3 for the direction of the axes. In this 
figure ¥ is the resultant stream velocity downstream of the shock 
(in the moving frame of reference), and the 4 and q axes are 
indicated. The final flow pattern depends cirucially on whether W is 
subsonic or supersonic; the criterion depends, in turn, on the Mach 
number corresponding to and on the wave inclination 0. 
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Boundary conditions . - The "boundary conditions Just downstream of 
the shock will now be obtained by application of the shock -wave relations. 

By geometry (fig. 4) the stream velocity components normal and tan- 
gential to the undisturbed shock are, respectively, 

cos e 


Y * sin 6 


The shear wave will provide directly a perturbation to and 

will cause indirectly a perturbation a(y) to the shock-wave angle, of 
initially undetermined magnitude. The effect of a is equivalent to 
an increment in 6. The associated perturbations to and V are 

found by obtaining their respective differentials and replacing d\fj^ 
by w^ and AO by a therein; the results are 

dU^ = cos 0 - aW^ sin 6 

dT = sin 6 + aWj^ cos 9 



The corresponding change in normal velocity U downstream of the shook 
is obtained from the normal-shock relation 


U 




1 + 


A TT 2 


By logarithmic 
temperature is 


differentiation and use of the fact that the upstream 


dUfl 


dUfl 


constant (whence — — b there is finally obtained 


\ 


U 


Ua 




( 4 ) 


where m = "Oa/u. 

On the downstream side of the shock the velocity perturbations in 
the directions of ? and q, respectively, are (fig. 4) 
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Wq = (U + dU) cos((p + ct) + (V + dV) sin((p + a) - ¥ 
= - (U + dU) sin (<p + c) + (T + dV) cos(q) + a) 


(5a) 


Equations (3) and (4) may "be used to evaluate the right-hand side 
of equation (5a), A first-order approximate result is obtained by 
taking cos a = 1, sin a a a and neglecting a tan (p and a cot cp 
in comparison with unity. It will be useful also to introduce the 
geometrical relation cos 9, the definition U/^/U » m, and 

to eliminate 6 by means of the oblique-shock relation 
tan <p = m tan 9. The final resirranged result is 



(5) 


These are the desired boundary conditions in a somewhat general form. 

In the present problem the perturbation w^ is associated with an 
incident sinusoidal shear wave parallel to (or to 

(figs. 2 and 3). It will be shown later that a refracted sinusoidal 
shear wave parallel to W (or to ^) will also arise. A suitable 
defining equation for w^ is 

— = e cos kriA (6) 

"A 

where k is the wave number (2rt/k *s wave length) . The coire spending 
etrgument for the refracted shear wave will Involve q and an altered 
wave number x. The argument of the upstream and downstream waves 
must match along the shock front, so that 

kT|A = along shock 

(by geometry (fig. 3), ^ = H P - g - iP . -nmg 
\ X cos 6 J 



e cos XT) 


along shock 


( 7 ) 
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Since the disturbance is sinusoidal^ the shock inclination a can like- 
wise be es^ected to be sinusoidal. For generality a phase shift can be 
allowed for, so that a can be assumed to have the fom 

a » e(a cos xt) + b sin xi]) (8) 

Substitution of these sinusoidal relations into the general form of 
the boundary conditions, equations (s), yields, after rearrangement. 


^ f 1 

eU [_m V 

- 2 m + m^'^ sin <P - f 1 - 2 2zi m\ cos <p + — — -■"1 

7+1 / \ 7+1 / cos (pj 

cos XT) + N 




'^o' r a 

eU [_■ m 

^ cos 1 ^ ^ “) ' 

^ COS XT) + 


[■ m 77l “) cos 1 ^ 

J 


(9) 


Equations (9) give, in final form, the conditions Imposed by the 
shock wave on the components peirallel to 4 and q, respectively, of 
the perturbation velocity immediately behind the shock; the parameters 
a and b therein governing the shock inclination a are undetermined. . 
These equations constitute the boundary conditions for the perturbation 
flow downstream of the shock. 


Evaluation of rotaticai term in governing equation . - Before 
equation (1) can be solved, the vorticity term (rotation term) on the 
right-hand side must be evaluated for the region behind the shock. A 
corresponding term has been evaluated in reference 5 for the flow 
behind a normal shock peirburbed by an isoenergetic upstream disturbance. 
This work has proved a useful guide, but it has been necessary to make 
modifications both for the variation in energy (that is, in total 
enthalpy H) and for the inclination of the shock in the moving 
frame of reference. The derivation is as follows: 


Downstream of the shock, the enthalpy H and the entropy s 

(and hence the vorticity) are constant along streamlines, and in the 

linear theory the streamlines are approximated by lines tj = constant. 

_ , 3s 

Thus, and gj. 


may be evaluated at the shock and the result will 
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hold downstream thereof if expressed as a function of t] alone 
( C eliminated ) . 


The total enthalpy upstream and at the shock is 
2 = Cp^A + i(WA + wa)^ 


“ Cp^A + 


I (> • ?) 


Hence, at the shock 


-S!- ^ -s; WLong Bhoct 


(10) 


The entropy upstream of the shock is constant by virtue of the 
assumption of constant pressure and density there. The entropy change 
in crossing the shock is given in terms of the upstream velocity hy 
(reference 6, eq.uation 144) : 


a 





(W^ + cos^(0 + a) 


- lip 
+ 1 “ 


(7 - l)(W^ + cos^(0 +a) + 

Cr + cos^(0 +<4 



Hence, on writing the differential and expanding the result under the 
assumption that is constant and and a are small, there 

is obtained 


6s 



(m 


1)2 (^ - a tan e) 


(11a) 


and 


y along shock 


^ T 





(11) 


Eecall now that the governing equation (l) reads 


(1 - 


-5i !!a 

W W 


where the right-hand side is the rotation term in question. The 
factors Hq and s^ have been evaluated in equations (10) and (11a), 

respectively; substitution with use of the geometrical relations 
of figure 3 yields 



10 


NACA TN 2864 


(1 - w2)+55 + 


Urn' 


2 cos cp 
cos^Q 



- U cos (p(in - 



( 12 ) 

where the right-hand side is to he evaluated along the shock (z = O) and 
expressed as a function of Tj alone. 

In the form (12) the governing equation has not yet been special- 
ized to a shear flow that is sinusoidal. The substitution of equa- 
tions (7) and (8) for respectively, introduces the 

sinusoidal character; furthermore, the relation tan (f> *= m tan 0 
can be used to eliminate 0; after simplification 

(1 . = Ue^-x Fseo <p + 2(m - 1) oos <p + a ^ sin sin xi) + xb ^ m ^ 


(13) 


Equation ^13) is the partial differential equation to be satisfied by 
the flow downstream of the shock subject to the boundary conditions 
( equations ( 9 ) ) . 


SOLUTION FOE HOEI^OHTAL WAVE 

The governing equation and boundary conditions have been set up 
for the general case of an inclined shear wave. It will be worthwhile 
to solve first, however, the much simpler special case of the horizontal 
shear wave. The results will illustrate important features of the 
general case as well as provide a limiting case of the general solu- 
tion, useful as a check. 

The horizontal wave is obtained by setting 0 » cp s 0 in the 
earlier equations; as a consequence T-»0, W -*-U, Ti->y, and 

x-»k. The governing equation reduces to 

+ tyy = - kUe(2m - l) sin ky (l4) 


vhere 


p2 = 1 _ ^ 


The boundary conditions "(equation (9)) reduce to 


2666 



2666 
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Ue 


Ue 




cos ky 


"^O 

= j]f~ = (m - l) a cos ky + (m - l) t sin ky 


12 

Ue 


(15) 


Particular Integral and complementary function . - A particular 
integral of equation (14) may be obtained ty inspection as 

tp = ^ (2m - 1) sin ky 

To ottain a complete solution there must he added a complementary 
function satisfying equation (14) with the right-hand side set equal 
to zero. The boundary conditions at x = 0 require that the function 
possess a sinusoidal variation with y. Such a solution will also 
contain an exponential factor, showing either amplification or atten- 
uation of the disturbance with distance x downstream of the shock; 
the case of amplification must be ruled out as physically unacceptable. 
These considerations limit the solution to the form 

- ^ 

= Uede sin ky 


where d is a constant of integration. 

The complete solution is the sum of \j/p and 






sin ky 


(16) 


Evaluation of undetermined constants . - The velocity components 
are obtained from equation ^16) as 


u = ty = Ue 


( -f) 

\2m - 1 + kde ^ j cos ky 


kx 


V = 


= pUekde ^ sin ky 


( 17 ) 


The undetermined constants and d are evaluated by setting 
X = 0 and comparing with the boundary conditions ^ equations (15)^ 
equating the respective coefficients of sin ky and cos ky« The 
results are 
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a =5 0 


■b 

d 


48m 

(7 + l)(m 
4m 

kCr + 1) 


ly 



(18) 


Velocity components . - Insertion of the value of d into eq.ua- 
tion (1?) yields the final result for the velocity components valid 
everywhere downstream of the shock 


u = Ue 


kx 

2m - 1 . e P 

7 + 1 


cos ky 


kx 

V = - pue — e P sin ky 
7 + 1 


all X = 0 (19) 


Just behind the shock 


'^o = (■! + 2 m^ cos ky 


> X = 0 


(19a) 


Vq = - PUe ^ 


and far downstream 


Uoo = Ue(2m - l) cos ky 


X = 


(19b) 


v„ = 0 


These results and the associated streamline pattern are exhibited 
pictorially in figure 5. 


These perturbation velocity components downstream of the shock 
are to be compared with the corresponding velocity components in the 
shear flow upstream of the shock (cf, eq.uation (?)) 


> 

= Uem cos ky _ 


Va = 0 
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The ratio of u»/u^ is 


— - 2 - 
- “ m 


( 20 ) 


Since m = order that a normal shock exist, it appears 

from eq.uation (20) that the normal shock always amplifies the 
horizontal shear wave, the maximum eunplification of ^ + 5 teing 

approached as the initial Mach number approaches infinity. 

Shock perturbation . - The local inclination of the shock from the 
vertical is given by equation (8) . With the previously determined 
values of a and b (equation (18)) inserted, eind with ky in place 
of xti, the inclination is 

48em ^ , 

(r V D fm - i) 


a = 


If the local shock displacement in the x-direction relative to the 
mean shock plane is called &x(y), then 


5x 




a dy 


4Pem , 

k(y + l)(m - 1) ^ 


( 21 ) 


Thus the shock, displacement curve is in phase with the velocity 
perturbation in the shear wave upstreeim of the shock (fig. 5). 


SOLUTION WHEN FLOW DOWNSTREAM OF EQUiVALENT OBLIQUE 
SHOCK IS SUBSONIC (W < l) 

The present case is a generalization from the horizontal wave 
Just discussed to a wave of eirbitrary inclination 6, The restriction 
to a subsonic mean velocity W behind the equivalent oblique shock 
insures a qualitative similarity of the flow: the governing equation 

is elliptic in both cases. Accordingly, the horizontal-wave result 
can serve as a guide. 

Governing equation and particular integral . - The governing 
differential equation (13) may be written in abbreviated form as 
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“ “ xU 6(A sin xt] - B cos xt) ) 


( 22 ) 


where 


cos (p + a ~ sin 
m 


A = ^sec (p + 2(m - l) 

B = h -fe.." sin <p 
m 

= 1 . 55 ® 

A particular integral is seen to he 


.] 


(23) 


(24) 


Complementary function . - From the result for the horizontal wave 
the complementary function should he expected to attenuate exponentially 
downstream of the shock, and from physical considerations the attenuation 
should depend upon the distance measured normal to the shock front, that 
is, upon X rather than, say, C* ®ie functional form that has the 
desired attenuation and possesses a sinusoidal behavior at the shock is 


'•'C 


cos <p - T 1 sin 
e 


<p)r sin 

Uos cos 


»)]} 


(25) 


where ? cos (p - il sin cp may he recognized as . just x. 


The arbitrary constant a in equation (25) is determined by a 
consideration of the boundary conditions (equation (9)): the argument 

of the cosine must reduce to xt] along the shock front, where 
? = T] tan <p. This requirement gives a = x cos (p/p^. Finally, when 
constants of integration c* and d’ are included, the complementary 
function is written as 


\|Tp = Uee 


XR 

— — cos q)(C cos <P - 11 sin (p) 

p2 




X cos <p 


|-2- (^ sin ip + T) COB + 

jji' sin - ^ (C sin <p + cos ip)j^ 


( 26 ) 
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Velocity components with undetemilned constants « - The complete 
solution for the perturbation stream function is 


t = 'kp + tc 


This expression (cf. equations (24) and (26)) contains four arbitrary 
parameters a and b (which occur in A and respectively, equa- 
tion (23)) and c' and d*, which remain to be determined. First the 
corresponding expressions for the velocity components will be obtained - 
they will be needed anyway - and then the boundary conditions on these 
velocities at the shock wave will be applied for the determination of 
a, b, c*, and d*. 


The perturbation velocity components in the direction of ^ and 
T) are w = and w* = - respectively; by differentiation of 

equations (24) and (26) there results 


- ^ 

^ = A cos XT] + B sin xq + p"^e 


cos tp(? cos <P - q sin <p) 


( 


X cos <p(C sin (p + q cos q>) 

(c sin (p + dp^ cos <p) cos r + 




(-cP^ cos (p + d sin cp) sin 


X cos <p(? sin cp + p^2 q cos <p) 




> (27) 


, “ — ^ cos <p(C cos <P - q sin <p) 

0e = P X 


( 


(cPvr^ cos (p - dPw sin cp) cos 


X cos (p(? sin {p + P^^ q cos <p) 


(cPw sin <p + dPw^ cos (p) sin 


X cos {p(^ sin <P + q cos <p) 


» 


where c' and d’ have been absorbed for convenience into new 
constants 


c = c’xPy. cos (p 


d = d’xp^ cos <p 
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The undetermined constants may now be considered as b, c, and d. 


Conditions along the shock on the downstream side have been 
designated by subscript zero; here % cos (p «* tj sin <P, euid the argu- 
ments of the exponential and sine and cosine terms reduce to zero 
and XT) respectively; 


^ = ( A + -% sin (p + Pw ’^’1 

Ue \ 32 p2 j 


[B - S- ^ cos (p + -% sin (p) 
V p2 ^ p2 J 


sin XT) 


(28) 


“W 

= P"2(cp,^2 gQ 3 <P _ cip^ sin (p) cos xt) + 


sin <p + dp^2 qos (p) sin xt) 


Evaluation of undetermined constants . - Equations (28) must 
agree identically with the boundary conditions (equations (9)) imposed 
by the shock wave on Wq and Wq*. Therefore the respective coef- 
ficients of sin XT) and cos xt) are to be equated; this yields four 
simultaneous equations for the four undetermined constants a, b, 
c, and d. In the reduction of the solution to final form certain 
alternative forms of the oblique-shock relations, given in appendix C, 
have been used. The results are 


a s= 


- CE + DF 
m ^ ^ 

C^ + d2 


m 


CF - DE 
C^ + 


c = - D* - F» 

m 


d = - D* 

m 


( 29 ) 


where 
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SOLUTION WHEN FLOW DOWNSTREAM OF EQUIVALENT OBLIQUE 
SHOCK IS SUPERSONIC (W > l) 

When the mean velocity W behind the equivalent oblique shock is 
supersonic, the solution must exhibit Mach waves. If the cross-stream 
velocity V of the moving reference frame is subtracted out, these 
waves appear to be moving downward (cross-stream) with the velocity V, 
If another transformation of axes is made so that the reference frame 
is "convected" downstream with the stream velocity U, then the Mach 
waves can be identified as plane sound waves moving normal to the wave 
fronts with sonic velocity. Mach waves and plane sound waves are, of 
course, the same phenomena viewed relative to different frames of 
reference . 


Gioveming equation and particular integral . - The governing 
equation (22) changes from elliptic to hyperbolic when W exceeds 
unity (that is, when W is supersonic). The particular integral is 
unchanged thereby and is still given by equation (24). It is found 
that the final solution yields b = 0 (and hence, B = O), and so it is 
convenient to delete the B-term at the outset; the particular integral 
is thus 


\|Tp 


UeA 


sin XT] 


Complementary function . - The complementary function satisfying 
equation (22) must be of the general form 

= f<C + p^q) + g(? - p^q) 


where 


p^= l\/w^ - 1. 
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The function f represents Mach waves inclined downward hy the Mach angle 
p from the ?-axis and the function g represents Mach waves inclined 
upward by the Mach angle. If attention is restricted to the range of 

shear-wave inclinations then the g-family of Mach waves can 

2 

be shown to represent disturbances overtaking the shock wave from 
behind. This property is related to the fact that, for a finite shock 
strength, the Mach angle is always greater than the angle between the 
shock and the |-axis. Since the disturbances actually originate at the 
shock wave by virtue of the passage therethrou^ of the initial shear 
wave, such Mach waves cannot arise, and the g-functlon must be zero. In 
what follows it will suffice to limit the discussion to the specified 
range 0-6— §, since the results for the remaining range 0 — 6 ^- J 

are readily obtained therefrom from symmetry considerations. 

The function f must reduce to 


f ~ sin XTi 


along the shock front, where 4 = il tan (p, in order to satisfy the 
boundary conditions (with b = O). A suitable complementary function 
is therefore 




Uec” , 

sin r: 

X 3w + ■tan <p 


where c*’ is a constant of integration. 


The complete solution for the perturbation stream function is thus 


I . . Ue 

\|r - tp + - — 


A sin xTi + c ’ ’ sin 


x(C + 

+ tan cp 


(31) 


This expression contains two arbitrary parameters a (occurring in A) 
and c’’ which remain to be determined. First the corresponding 
expressions for the velocity components will be obtained, and then the 
boundary conditions on these velocities at the shock wave will be 
applied for the determination of a and c ’ ’ . 


Velocity components with undetermined constants , - The perturba- 
tion velocity components in the direction of ^ and t| are w = f 
and w* = respectively; the expressions are 


w 

U6 


A 


cos XT] + 


c sec (p 
P^ + tan (p 


cos 


x(C + 

P^ + tan (p 


p^c sec (p 
^ “ P^ + tan (p 


P^ + tan <p 




( 32 ) 


COS 
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where the constant c’* has heen ah sorted into a new constant 
c = c • * Py cos cp . The undetermined constants are now a and c . 


Along the shock C = "tan (p, and the arguments of all cosine 
terms reduce to the expressions for the velocity components w 

and w* he come 


^o + c sec cp \ 

Ue = + tan cpj 


cos XT) 





P^ c sec cp 

•T 7 cos XT) 

p,^ + tan (P 


(33) 


Evaluation of undetermined constants , - Equations (33) must agree 
identically with the "boundary conditions (equation (9)) imposed hy the 
shock wave on w^ and Wq*» If the respective coefficients of sin xr) 

and cos xt] are equated, there results h = 0 and two simultaneous 
equations for a and c. Thus, the initial specification of h = 0 
has heen Justified a posteriori. 


The solutions may he written in the form 

C* + GF‘ 


a = m 


E‘ + GD’ 


c = - D* - F‘ 
m 


(34) 


where 


c = 

2 21 ^ m-2ri + (m- l) 


r + 1 

L 

D* = 

(m - 1) [ 

1 + (m - l) cos2(^ 

E' = 

(m - 1)2 

sin (p cos cp - (l + 

F' = 

2(m - 1) 

sin cp cos cp 


1 - P tan cp 

G = —7 — = tan (|i - cp) 

P^ + tan cp ' 


(35) 


where pi = cof^P^^ is the Mach angle. (The definitions for D* and 
F’ herein are unchanged from those included in equation (30).) 
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RESULTS ARD DISCUSSION 
Velocity Field 

Ttie Telocity field downstream of the shock wave, produced hy con- 
vection of an ohliq.ue sinusoidal shear wave through the shock, has been 
calculated; the results are distributed throu^ the preceding sections. 
The main results will now be presented in more compact form,, simplified 
to aid in the geometrical interpretation. (The special case of the 
horizontal shear wave was discussed earlier.) 

Frames of reference . - The analysis has been carried out in a 
special frame of reference in which the flow is steady; all formulas 
will be given relative 'to this steady-flow frame. Also of considerable 
Interest is a frame of reference convected by the mean flow downstream 
of the shock; this frame is at rest relative to the general mass of 
fluid there. The relation between the two frames is shown in figure 6. 
Fonnulas relative to the steady-flow frame may be converted to apply 
to the convected frame by means of the transformations 

+ Wt 

+ Ut 
y-*-y + Vt 

The criterion on. W . - Although the stream velocity U downstream 
of the specified normal shock (fig. l) is always subsonic, the nature 
of the flow depends primarily on the stream velocity W downstream of 
the equivalent oblique shock (figs. 2 and 3), which may be either sub- 
sonic or supersonic. The velocity W may also be interpreted as the 
relative velocity of the steady-flow frame of reference and the con- 
vected frame (fig. 6). Two forms of the solution for all flow quan- 
tities thus appear one for the subsonic ran^^e W < 1, the other for the 
supersonic range W >1. The dividing line W = 1 is what has been 
designated "the criterion on W" at th£ head of this section. Since 
W depends on the initial Mach number and the inclination 0, the 

equation W = 1 g^ves, in effect, a relation between a critical 
value of 0 and U^. The relation is conveniently expressed in terms 

of m = U/^A^, which depends on (see appendix C) : 

0 ( 31 * — i tan ^ 
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A graph of j®cr| ■'^®rsus 

Resultant velocity, W<1 « 


is given in figure 7. 

- Eq.uations (27) may he recast in the form 


= S cos jjty(y - X tan (p) + 6^ + n(z) cos ^Xy(y - z tan <p*) + 8^ 
= Pwn(x) sin |»y(y - I tan <i)‘) + 8pJ 


(38a) 


where 


jw^l = W^e = amplitude of sinusoidal velocity vp^ in initial shear wave 
Xy = X cos <p = k cos 6 

g = co^ ^|a2 + b2; a - A(a), B = B(h) 

cos e "^yPw/P^ 

n(z) H ^ e 

5g = tan-1 (^) 


<p« = 



- tan 
tan“l 


_1 tan <P 

p2 

cP^ - d tan cp 
dp^ - c t£in (p 


The functions A and B . are given hy equations (23) and a, h, c, and 
d are the Initially undetermined constants which have heen evaluated in 
equations (29) and (30). 

Resultant velocity, W>1 « - Equations (32) may he recast as well 
as generalized to apply for both positive and negative values of (p 
as follows: 


, ^ -r = S COS X (y - X tan (p) + n cos x (y - x tan (p*) 
|Wa| y ^ 

Jil— = P^ n COB X (y - X tan <p‘) 

I^a| "" y 


V (38b) 



c 
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where 

Iw^l H W^e = amplitude of sinusoidal velocity w^ in 
' initial shear wave 

s E £22_« A; A . A(a) 
m 

Xy = X cos <P = k cos 0 

_ _ cos 0 sin u 
~ m sin(ii + (p) 

<p* = angle of magnitude (|<p| - and has the same sign as cp 
p = Mach angle = cot”^P^ 

The function A is still given hy eq.uations (23), and a and c are 
evaluated in equations (34) and (35) o 

Shear-wave component . - The cosine in the S-term is constant along 
lines y - X cot cp = constant; such lines are inclined at an angle cp 
with the horizontal and are thus parallel to the ^-axis. Since w 
is parallel to 5 and- w* is parallel to ti, it is seen that the 
S-term represents a pure shear flow parallel to the C~axis. Stated 
otherwise, this is a rotational flow; the rotation (or vorticity) is 
just 2 , which was evaluated earlier in terms of gradients of entropy 
and total enthalpy (cf. equations (l) and (13)). The shear flow may 
be described also as an incompressible, plane, transverse, sinusoidal 
wave. 

The amplitude and phase of the shear wave are compared with those 
of the initial shear wave in figure 8 for an initial Mach number of 1.5. 
The amplitude amplification ratio is S and the angle of phase lead 
is 5g; both are plotted against the initial wave inclination 0. There 

is seen to be a small phase lead in the_ subsonic range (W<l) and 
none at all in the supersonic range (W>l). The amplification is 
nowhere_less than unity, with a cusp-like peak of 1.73 at the sonic 
point ¥ = lo 

Pressure-wave component . - The remaining terms in equations (38a) 
and (38b), involving the factor 11, correspond to an irrotational 
velocity field, or potential flow. That is, if the derivation is 
traced backward, the II -terms are found to have come from the complemen- 
tary function, which is a solution of the governing equation with the 
vorticity 2 set equal to zero. This part of equations (38a) and 
(38b) defines what may be called a pressure wave since there is 
associated with it a first-order pressure field : the shear wave 

contributes nothing to the pressure. 
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The pressure wave may he interpreted as a distribution of sound 
waves. This interpretation is particularly evident for the case 
W>1^ where the solution has been obtained in the form of Mach waves;, 
if a transformation is made from the present special frame of refer- 
ence, relative to which the flow is steady in time, to a frame moving 
with the general stream, then the Mach waves will reappear as plane 
sound waves moving normal to themselves with sonic speed. 

The same transformation results in somewhat more complication 
when W<1: the resultant pressure pattern does not then propagate 

with the speed of sound, but it can be represented (as can any 
two-dimensional irrotational gas-flow field) as a superposition of 
cylindrical sound waves which individually propagate with sonic speed. 
The associated velocity pattern in this case exhibits the following 
features, which are brought out by an examination of eq.uations (38a): 

The radius vector in a graph of w’ versus w (hodograph) moves in 
an ellipse when x is held fixed and y varied; the major and 
minor axes are n|w^| and respectively. At x = 0 the phase 

angle relative to the incident shear wave is 6^. On the other hand 

the argument of the cosine and sine is constant along lines 
y - X tan <P* = constant; these are lines inclined at an angle cp’ to 
the horizontal. Along such lines the perturbation velocity (w,w*) 
remains constant in direction but attenuates exponentially with x; 

the exponent is 

p2 

For the case W>1, the velocity pattern associated with the 
pressure wave is much simpler (equations (38b)). The perturbation 
velocity vector (w,w') is constant along lines y - x tan <p* = 
constant and is, in fact, normal to such lines. In this case 
(p' = (cp - Mach angle), and these are just the Mach lines (or 
envelopes of the sound waves); they are inclined downward by the Mach 
angle p relative to the C-axis. It will be noted that the 
definitions of <p'j_ the inclination angle of the lines of constant 
phase, a^ee at W = 1, although expressed differently for W <1 
and for W > 1 . 

The amplitude and phase of the w and w’ components of the 
velocity in the pressure wave are compared with the amplitude of the 
initial shear wave in figure 9 for an initial Mach number of 1.5. The 
amplitude amplification ratios are H and respectively; FI, 

pi^n, and a phase angle (lead) Sp are plotted in the curves against 

the initial wave inclination 6. In the subsonic range (W<l) FI 
and P^Fl attenuate exponentially with x and only the values for 

X = 0 are plotted. The phase lead varies from 180° to zero in this 
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subsonic range and remains zero throughout the supersonic range (W>l), 
A rather striking feature is the relatively small perturbation velocity 
in the supersonic range. Thus, although the incident shear wave can 
give rise to a simple sound wave upon passing into the shock wave, the 
particle velocity in this sound wave amounts for most cases to 
2.0 percent or less of the velocity in the initial shear wave for 
% = 1.5. 


Pressure Field 


It is shown in appendix B tliat the perturbation pressure is 
related to the velocity according to eq.uation (Bll); in the present 
notation this becomes 


or 


8p = - pWWp 




(B11‘) 


Here Wp is that component of the perturbation velocity associated 

with the pressure wave and directed parallel to W (that is, along the 
5 -axis). Equation (Bll*) may be recognized as the linearized Bernoulli 
equation as limited to the velocity in the pressure wave. 

Upon substituting for W and W and using for Wp equations (38a) 
and (38b) with the S-terms omitted, there results 

^ ‘ <7S)m [’V<y - ^ * ®p] 

where 6p is to be taken as zero in the supersonic range of ¥. 

This result for the perturbation pressure is proportional to 
n sec (p; n has been plotted in figure 9, together with 6p, as a 

function of wave inclination 0 for = 1.5. 


Density Field 

The density perturbation is related to the velocity and entropy 
perturbations according to equation (B12) of appendix Bj in the 
present notation this is 
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P 


W ■ Cp 


(B12‘) 


The term in Wp is the contrihutlon of the pressure wave. This term. 

differs from 5p/p (eq.uation ^7 a simple factor I/ 7 , so that 

the contrihution is obtained at once from equation (39). 


The term in 8e is the contribution of the shear wave. The 
entropy perturbation &s has not been given explicitly before, but it 
can be obtained from equation (lla) by use of geometrical relations 
and the known result for a (see following section). Upon evaluation, 
the term in 8s is found to be 


8s _ I^A. I 2(m - 1 )^ cos $ 

Cp Ua 2-±A m - 1 
7-1 


^a tein 0 - l) cos xti + b tan 6 sin xqj 


(40) 


Shock-Wave Perturbation 

The local perturbation in the shock inclination angle may be 
written (cf. equation (8)) 

a = e(a cos Xyy + b sin Xyy) 

where a and b are evaluated in_equation (29) for W<1 and equa- 
tion ( 34 ) for W>1 (b = 0 for W>l). 

The local shock deflection 8x from the plane 2=0 is obtained 
by Integration of the elope a : 



The result may be put in the form 


= ^aI °os(V 

where ^shock “ tan"^ is the phase angle and X = 2jt/k is the 

wave length of the initial shear wave. 

For a given wave length the factor ^/a^ + b^ is proportional to 
the amplitude of this sinusoidal corrugation in the shock wave; 
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Va^ + Is plotted against the initial wave inclination 6 in 
figure 10. The phase angle ^ahock also plotted: the shock- 

wave corrugation is in phase with the initial shear wave (^shock “ 

when the initial wave is horizontal (6 «= O) , The shock corrugation 
progressively lags the initial shear wave as 6 is increased until 
the sonic condition W *= 1 is reached; at this point the lag is 90°, 
and this value is maintained throughout the range W>1 as 6 is 
increased to 90° (v ertical initial shear wave). At 6 b 90° the 

amplitude factor Va^ + has fallen to zero: a vertical sinusoidal 

shear wave passing hy convection through a vertical shock wave causes 
no perturbation of the shock form or position. 


Intensity of Sound Field 

The acoustic intensity of the noise or sound field generated by 
the Interaction of the shock wave and the turbulence has been found 
to be relatively hi^.l It will suffice for an order-of -magnitude 
estimate to replace the turbulent field by a single plane wave, or 
Fourier component, with the same kinetic -energy density. Hourly 
this Implies that the root -mean- square turbulent velocity is to be 
identified with 0.707 j'WA|« 

The sound pressure is proportional to n sec (p, where n is 
plotted in figure 9. The relatively high values indicated for the 
subsonic range attenuate rapidly with distance x downstream of the 
shock; when x appreciably exceeds several wave lengths the values 
are negligible compared with those in the supersonic range. A rough 
average over all wave incli nations 0, assuming the subsonic range 
contributes nothing, gives |n sec (p| = 0 . 082 ; this value will be 
used in the noise estimate. 

The noise level in decibels relative to the standard reference 
level 6p^ = 2.015X10"^® atmospheres at standard density and speed 

of sound Sq is given by 


^A high acoustic intensity does not imply a large fractional 
pressure perturbation in relation to ambient pressure: the sensi- 

tivity of the ear is so great that extremely small pressure pertur- 
bations correspond to very loud noises. 
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db = 20 


= 20 

where the 5p’s are root -mean- square values. By equation (39) 

^ Vf ym . In sec <p I (42) 

p % (y + l)m - (r - 1) 

As an example the noise level generated by the turbulence passing 
through a normal shock in a representative supersonic wind tunnel will 
be estimated. A root -mean- square velocity of turbulence of 1 foot per 
second is assumed to exist in the test section where the mean speed is 
1400 feet per second and the Mach number is 1.5 = 1.5). 

Thus 0,707 Iwft^l and % are taken to be 1 and 1400 feet per second, 
respectively. A summary of these and the remaining parameters of the 
example is 

I n sec <p I = 0.082 
0.707 |w^ 1 = 1 foot per second 
= 1400 feet per second 
m = 1.862 (- % = 1.5) 

7 s5 1.4 

p = 0,272 atmosphere ( ~1 atm. reservoir pressure) 

6Pq c 2.015X10"^*^ atmospheres 


pa 

The estimate based on equation (42) gives a pressure perturbation 
5p rms/p = 7.50 x 10"5, and by equation (41) the corresponding sound 
intensity is 102 decibels. This represents very intense noise, reaching 
a level which can damage the ear on continued exposure (reference 7). 
This noise estimate is thought to be conservative, corresponding to a 
supersonic wind tunnel with a relatively low level of turbulence. It 
appears probable that many tunnels will considerably exceed this level. 


log + 10 log 

if 4 ) * ^ 
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Generalization to Oblique Shocks 

The analysis refers to flov throu^ a normal shock, but the 
results are easily generalized for oblique shocks. In the oblique- 
shock case the component of the upstream velocity normal to the shock 
plays the role of the component parallel to the shock is iguored 

in formulating the equivalent steady-flow problem. A formal approach 
is to retain the present definitions wherein TJp^ is the actual upstream 

velocity (taken horizontal) and 6 and <p are referred to the 
horizontal; the oblique shock is assumed inclined by some angle a 
measured clockwise from the vertical. Then the present formulas 
will be generalized to apply to the oblique shock if the following 
transformations are made; 


cos a 

e-*‘9 + a 
(p -*■ <p + a 


Related Problems 

The sound field produced downstream by the convection of turbulence 
throu^ a shock has been discussed. Also of Interest are sound fields 
incident upon a shock in the absence of turbulence. The elementary 
sound disturbance is the plane sinusoidal wave; a longitudinal wave. 

The passage of such a wave through a shock, which is an unsteady- 
flow problem, can again be converted to an equivalent steady-flow 
problem by transformation to a reference frame moving with a suitable 
velocity parallel to the shock front; in this frame the sound-wave 
pattern will appear as a stationary Mach wave pattern. A diagrammatic 
construction is shewn in figure 11, Note that either of two sound 
patterns of uniquely related incllnatians may be rendered stationary 
by a given choice of V; the two patterns may be identified with the 
two families of Mach waves in a stream of supersonic velocity W^. 

The equations for the boundary conditions at the shock and the 
vorticity behind the shock will be modified from those for the 
present case of the shear wave, but the general character of the 
solution will be unchanged. Thus, a shear wave as well as a sound 
wave will appear downstream of the shock. The discussion will be 
carried no fuirther here; the solution has been obtained in the 
unpublished investigation of F, E. Moore by hie unsteady-flow method. 

The interaction of a sinusoidal Mach wave pattern with a normal 
shock constitutes a simple special case; here the velocity V of the 
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moving reference frame may "be taken to te zero. This problem has teen 
solved in general terms ty Adams (reference 5); he limited his 
discussion, however, to the vicinity immediately downstream of the 
shock. The character of the flow further downstream can he inferred 
from the parallel that exists between this problem and the problem 
herein of the horizontal shear wave: in both cases Y is zero. The 

asymptotic flow far downstream is therefore a horizontal sinusoidal 
shear wave. Near the shock the wave is modified by transverse and 
axial components (with associated pressure perturbations) which attenuate 
exponentially with distance downstream of the shock (cf. fig. 5). 

According to these considerations, sinusoidal corrugations in a 
wind-tunnel wall, or a plate, upstream of a plane shock wave will 
generate a horizontal sinusoidal shear flow. Such a shear flow mi^t 
have applications in special experimental work. 


CONCinDIKG EEMARBS 

The effects produced by the convection of an Inclined plane 
sinusoidal shear wave through a normal shock have been analyzed. Such 
a wave may be interpreted as a single spectrum component of a turbulent 
field; that is, the turbulent field can be represented as a super- 
position of such shear waves of all orientations and wave lengths 
(Fourier Integral) . 

When the turbulence is convected through a shock, the Individual 
waves do not mutually Interfere if, as specified herein, the Intensity 
is sufficiently low; thus the modified field downstream of the shock 
can be obtained in principle by superposition of the modified indi- 
vidual waves. In practice the initial wave distribution is known only 
statistically, and statistical changes only can be calculated. In 
either case the present analysis of the behavior of a representative 
individual wave constitutes a prerequisite to the determination of the 
changes in the weak turbulent field. 

It is found that a sinusoidal shear wave of arbitrary inclination 
as it passes into the shock gives rise downstream to a shear wave of 
altered inclination and altered amplitude. In addition, there is 
generated a "pressure wave": an additional velocity field with 

associated pressure disturbances that can be recognized as sound waves. 

The analysis is made in a frame of reference moving with a certain 
velocity W referred to axes at rest relative to the general mass of 
fluid downstream of the shock; W is the vector sum of the reversed 
downstream velocity and the cross-stream speed of the ripple pattern 
in the shock wave. The results depend crucially on whether W is 
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subsonic or supersonic; wben W is subsonic both the shear wave and 
pressure wave are shifted in phase relative to the Initial shear wave, 
and the pressure wave shows an exponential attenuation downstream, of 
the shock; when ¥ is supersonic there are no phase shifts, and the 
pressure wave takes the form of a plane, undamped, sinusoidal sound 
wave. 


A weak initial shear wave is found to produce a surprisingly 
intense pressure wave or sound field downstream of the shock, as 
measured in acoustic terms. This implies that the convection of 
relatively low- intensity turbulence throu^ a shock will generate a 
very intense noise field in the downstream region. In an example 
the noise level generated by turbulence in a representative supersonic 
wind t unn el was estimated to be of the order of 100 decibels. 

Lewie Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, September 26, 1952 



APPEHDIX A 


SIMBOLS 

function defined in equation (23) 

parameter in shock-vave perturbation (equation (8)); also 
speed of sound 

function defined in equation (23) 

parameter in shock-wave perturbation (equation (8)) 

function defined in equation (30) 

function defined in equation (35) 


parameter 


/= c‘*8^^ cos q> for W>l\ 

\= cos <p for W<1/ 


constant of integration 
constant of integration 
specific heat at constant pressure 


function defined in equation (30) 
function defined in equation (35) 
parameter (= cos ff) 

constant of integration 
function defined in equation (30) 
function defined in equatiai (35) 
function defined in equation (30) 
function defined in equation (35) 
function defined in equation (35) 
stagnation enthalpy (per unit mass) 

wave number of shear wave in region A (incident shear wave 
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M Mach number (u/a, appendix B) 

m velocity ratio across normal shock (Vp^/u) 


P 

S 


s 

T 

t 

JJ,Y 

U 

u,v 

W 

W 


^,7 
P = 

Pv,= 



pressure 

relative amplitude of refracted shear wave (see equa- 
tions (38a) and (38b)) 

entropy (per unit mass) 

temperature (absolute) 

time 

stream velocity components in x- eind y-directions (fig. 3) 
(equivalent steady-flow problem) 

Mach number associated with U (u/a) 

perturbation velocity components in x- and y-directions, 
respectively, (fig. 3) 

stream, velocity in 4 -direction (resultant of TJ and V) 
(equivalent steady-flow problem) 

Mach number associated with ¥ (w/a) 

perturbation velocity components in?- and ij-directions, 
respectively (fig. 3) 

that part of w associated with pressure wave 
rectangular coordinates (fig. 3) 



ratio of specific heats 

phase lead of pressure wave relative to incident shear wave 

phase lead of refracted shear wave relative to incident 
shear wave 
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e 

e 

^cr 

X 

n 

p 

a 

<P 

cp' 

t 

% 


measure of strength of Incident shear wave (w^/w^) 

inclination of lines of constant phase in incident shear wave 
(figs, 1, 2j Z, and 4) 

critical value of 6 for which ¥ « 1 (function of U^) 

wave number of refracted shear wave 

Mach angle associated with W (sin”^(l/¥)) 

inclined rectangular coordinates (fig. 3) 

relative amplitude of velocity component v in pressure 
wave (see equations (38a) and (38b)) 

fluid density 

perturbation in local shock angle (fig. 4) 

inclination of lines of consteint phase in refracted shear wave 
inclination of lines of constant phase in pressure wave 
perturbation stream function 
complementary function (component of \|r) 

particular integral (component of \|r) 
vorticity (v^^ - Uy) 


Subscripts; 

A region A (upstream of shock) 

o evaluated at shock, on downstream side 

7} indicate the corresponding partial derivatives (e.g., 
n V 




; an exception is Xy = x cos <p = k cos 0 


(Unsubscripted velocity components, pressure, and density 
refer to region downstream of shock.) 


Prefix 

6( ) increment in ( ) 
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APPENDIX B 

LINIAR PERTUEBATION THIOEY FOR ROTATIONAL FLOW 

The generalized governing equation for the stream function can he 
obtained by extending Sears* consteint-energy development (reference 4) 
to include the effects of variation of energy (total enthalpy) o A 
different approach is, however, employed herein. Equations for the 
pressure and density fields are also obtained. 

In applying the results of this appendix to the developments in 
the main text it is to be noted that the x- and y-axes herein will go 
over, respectively, into the K- and q-axes therein; this is a conse- 
quence of the difference in direction of the main stream in the two 
cases. There is a corresponding change in the notation for the 
velocity components. 


Governing Equation 


Basic equations . - Consider the steady two-dimensional adiabatic 
flow of an inviscid fluid with local velocity u*,v*, pressure p, 
density p, temperature T, and entropy s. Assume only small pertur- 
bations from a uniform horizontal flow such that u* = U + u, v* = v, 
with u/c, v/lJ<<l, and also 5p/p, 8p/p, etc. «1. Then the basic 
flow equations may be linearized by neglecting quantities of order 
u/U, and so forth, in comparison with unity. A convenient form of 
these linearized equations is 


Continuity: 


State : 

P 


Energy: 

Momentum: 


+ v^ 


1 Dp 
P Dt “ 


Sp 5s 


pa 


Dt 


0 


-Px = PUUx 1 

-Py = I 

Ml 


0 


(Bl) 

(B2) 

(B3) 

(B4) 


where D/Dt signifies the Lagrangian operator for differentiation 
following the fluid motion. 


Elimination of density from continuity equation . - The Lagrangian 
form of the state equation is, by virtue of the energy equation. 
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1 DP D£ 

P Dt - p^2 I>t 

' [(U . U)p, . 

Upon linearizing, assuming Py and p^ to "be of comparaUle magnitude, 
this is 

1 Dp 

P Dt “ pa2- 

and by use of the first momentum equation 


i - - M^u 

p Dt ■ ” 


The linearized continuity equation (Bl) may 

(l - M^)u2- + v„ = 


accordingly be written 
0 (B5) 


Formulation of governing equation , - Define a stream function \|r 
such that 


U = ty 

V = - (1 - 


(B6) 


Then equation (B5) is identically satisfied by u and v as defined 
in equation (B6). The governing equation for \|^ is now obtained by 
expressing the vorticity v^j- - Uy = 2 in terms of : 

(1 - + tyy = - 2 (B7) 

A useful expression of the vorticity in terms of gradients of 
entropy and total enthalpy is given in reference 8, equation (8.3), as 


where q is the resultant velocity and h/bn. indicates differentiation 
normal to a streamline. In the small-perturbation flow the streamlines 
are approximated by the lines y = constant, so that b/b(\ = b/bj; also 
q = U. Thus 
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The governing equation for equation (B7), can now he amplified 
to read 

(1 - - 2 » (B8) 


This equation and its companion 


U » MTy 

V = - (1 - 


(B6) 


constitute the simplified generalization of Sears ' governing equation 
for linearized rotational flow (reference 4, equations (12) and (15)); 
Sears* equation is restricted to flows of constant total enthalpy H. 

Equation (B8) exhibits the following very interesting property: 

In the small-perturhation velocity field considered here the effect of 
the rotation or vorticity 2 is independent of how it arises, whether 
from a gradient of entropy or a gradient of total enthalpy, or a 
comhlnation of both. And it is only through their ■ contribution to a 
(and perhaps to the boundary conditions) that variations in entropy 
and total enthalpy affect the velocity field at all. ' 


Pressure Field 


Equations relating the pressure distribution to the velocity 
distribution will now be derived: The momentum equations (B4) may be 

rewritten in the form: 


since 


2 = Vj - Uy. 


- + Uu^ = 0 
^ + UUy » - U2 



(B9) 


Wow let consideration be limited to special types of flow such that 

ussu’+u*’"l 


V = V* 


> 


(BIO) 


KACA TN 2864 


37 


where u%v* is an Irrotational flow (v’ - u* = O) and 

iX T 

u" = u*‘(y) is a pure shear flow parallel to z (v” = O) . Then the 
vorticity 2 is given hy 

- Uy = - U“y 

and 

Thus equations (B9) hecome 


— + Uu»2: = 0 
P ^ 


h. 


+ Uu»^ = 0 


These two equations are equivalent to 

^ + Uu» = 0 


(Bll) 


which is just the linearized Bernoulli equation in terms u* alone. 

The physical interpretation is this ; If the assumed total pertur- 
hatio^consists of a plane she^ flow (u**iO) and a potential flow 
(u*,v*), then there is, no pressure perturhation associated with the 
shear flow; the entire pressure perturhatlon arises from the potential 
flow and is related to u* hy the ordinary linearized Bernoulli 
equation. In other words, the pressure is obtained hy subtracting out 
the shear-flow velocity and applying the linearized Bernoulli equation 
to the remaining velocity . 


Density Field 

The density distribution can be related to the velocity and entropy 
distributions as follows: The starting point is the differential 

equation of state (B2) 


hp _ 8p Ss 


pa^ 


Since H and s are constant along streamlines (reference 8), this 
approximates the general small -perturbation flow to the extent that the 
lines y = constant approximate streamlines. 
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Again assume that the flow is a combined potential flow (u*,v’) and 
shear flow (u'',0).^ (See equation (BIO) and after.) Then equation (Bll) 
applies for 5p/p^ and the density field is given by 


or 


= _ Hlil _ ^ 

P a^ °p 


6P m 2 ]il ^ 
P " U ■ Cp 


(B12) 


Thus it is found that the density perturbation depends on the potential 
flow via the velocity perturbation u* and on the shear flow via the 
entropy perturbation 5s . 


Since H eind s are constant along streamlines (reference 8), this 
approximates the general small -perturbation flow to the extent that the 
lines y = constant approximate streamlines. 
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APPMDIX C 

VARIANTS OF THE SHOCK EEIATIOIK 

The ratio of the normal velocities before and after the shock has 
been defined as m 


m , U^/U 

Thus by reference 6, eq^uation 0-14), 


(Ql) 


m = 


y + 1 n 2 
2 ^A 


1 + 


n 2 


(C2) 


where is the normal Mach number ahead of the shock. Correspondingly, 

U is the normal Mach number behind the shock, and by reference 6, 
equation (112), 


1 -f y ~ ^ U 2 

ij2 = 

yU 2 _ 7 1 . 

^A 2 


From equations (C2) and (C3) it can be shown that 


(C3) 




j2 (y + l)(m - 1) 


(C4) 


and 


V - 1 


1 * ^ 


= m - 1 


(C5) 


where 


3 


^ = 1 - TJ^ 


The equality of transverse velocity components across an oblique- 
shock wave requires, in the present notation, that 


tan e = U tan cp 


Then, with the definition (Cl), 
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m tan 9 = tan <P (C6) 

Equations (C^) and (C6) together allow (p to be determined in terms 
of 6 and U^. 
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Shock 

Figure 1. - Convection of plane oblique sinusoidal shear wave 
through shock: original unsteady -flow problem. 



Figure 2. - Transformation to equivalent steady-flow problem by 
superposition of velocity V. 
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Shock 

Figure 5. - Passage of horizontal shear wave through normal shock, showing 
perturbation of shock and final amplification of shear wave. 



Figure 6. - Relative motion of reference frame moving with general 
downstream flow (convected frame) and reference frame of analysis 
(steady-flow frame). The steady-flow frame moves downward along 
the shock front with a component velocity V and carries the 
ripple pattern with it. 
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Figure 7, - Upper curve shove variation with initial Mach numher of critical wave inclination for 
which W is sonic. .Lower curve shows variation of m = U^/U “ tan (p /tan G with Initial 
Mach numher. 
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Figure 10. - Amplitude and phase of ripples developed in shock by passage of shear wave. 

Initial Mach number 1.5. 
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Figure 11. - Construction for translation V to render either of two 
sound-wave patterns stationary in a main stream U^. 
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